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Abstract. We studied the formation and evolution of low-mass stars within halos with high concentration of dark matter 
(DM) particles, using a highly sophisticated expression to calculate the rate at which DM particles are captured inside the 
star. For very high DM densities in the host halo (p% > 10 10 GeV cm~ 3 for a 1 M,?, star), we found that young stars stop 
sooner their gravitational collapse in the pre-Main Sequence phase, reaching states of equilibrium in which DM annihilation 
is their only source of energy. The lower effective temperature of these stars, which depends on the properties of the DM 
particles and DM halo, may be used as an alternative method to investigate the nature of DM. 
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Cosmological observations at different scales suggest that ~ 96% of the Universe is composed of invisible matter and 
energy. From that ammount, ~ 23% is known as Dark Matter (DM), which is believed to be partially responsible for 
the formation of the first structures of the Universe. If the DM component is composed by particles, they are very 
likely WIMPs (Weakly Interacting Massive Particles), given that they must be stable, neutral, and massive. WIMPs 
have sizeable scattering cross sections with baryons, which upper limits are set by direct detection experiments [1,2], 
and their annihilation cross section is known to be of the order < o a v >= 3 • 10~ 26 cm 3 s _1 [3]. These two properties 
make DM particles accumulate inside stars and annihilate among themselves, providing a new source of energy for 
the star [4, 5]. 

Recently, many authors have been studying the influence of self-annihilating DM particles in the first generation 
of stars [6, 7, 8, 9, 10, 11, 12], in compact objects [13, 14, 15], in the Sun [16, 17, 18, 19], and in low-mass stars 
of the Main Sequence (MS) [20, 21, 22]. In this proceedings we will review the new stellar evolution scenarios that 
low-mass stars follow when formed in dense halos of DM, which are described in detail in our work [23], but this time 
we will use an upgraded calculation of the process of capture of DM particles by the star which, compared with our 
previous results, leads to stronger effects in the forming star in the pre-MS phase. Nevertheless, the overall results and 
conclusions are similar. 



The evolution of stars in DM halos is regulated by the capture of DM particles in the stellar interior. DM particles 
in the halo may scatter with the nuclei of elements inside the star. If in these collisions the DM particles lose enough 
energy, they will get gravitationally trapped in the core of the star. To compute the number of particles that are captured 
per second by the star, we adapted part of the publicly available DarkSUSY code [24], which calculates the capture 
rate C x from the original expressions of A.Gould [25]: 



where Q. v (w) is the probability of a DM particle with a velocity w to have, after the collision with the nucleus of an 
element i, a velocity v lower than the escape velocity of the star v esc ,r at the radius of the collision, 
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TABLE 1. Capture rates obtained for stars at different stages of their evolution, 
for different stellar masses, stellar velocities, and mass of the DM particles. If not 
stated otherwise, the capture rates are those when stars are in the Zero Age Main 
Sequence (ZAMS), with Af* = 1 M Q , v* = 220 km s _1 , and m x = 100 GeV. In 
all cases, the stars evolved in a halo of DM particles with p x = 10 8 GeV cm~ 3 , 
v x — 270 km s _1 , and (J x .sd = 10~ 38 cm 2 . 



C x (s- 1 ) 



Stage 


Pre-MS* 


1 x 10 31 


ZAMS 


9 x 10 32 


M* (M s ) 


0.7 
3.0 


4 x 10 32 
8 x 10 33 


v* (kms -1 ) 


50 


2 x 10 33 


350 


3 x 10 32 


m x (GeV) 


10 
1000 


3 x 10 34 
1 x 10 31 



"(when L = 10 L 0> T eff = 4610 K) 



„_m x _/i±l m 

m n 2 

In our computations we calculated the capture rate for 16 elements (H, 4 He, 3 He, 12 C, 14 N, le O, Ne, Na, Mg, Al, 
Si, S, Ar, Ca, Fe and Ni), which abundances are followed by our code. The initial chemical abundances and initial 
metallicity of the stars were assumed to be as the solar ones [26]. All elements except hydrogen have negligible 
contributions to the total capture rate. This is because H is the only element with spin dependent (SD) interaction with 
the DM particles, and the experimental scattering cross section limits for this type of interaction is less stringent than 
for the spin independent (SI) one. 

The calculation of the capture rate implemented in the present work takes into account that the DM particles have a 
velocity distribution fo(u), which we assumed to be a Maxwell-Boltzmann distribution with a dispersion velocity v~ x . 
This C x expression also permits the consideration of different velocities of the star v*, leading to a distribution of the 
velocities of the DM particles seen by the star / Vjr (u) equal to: 



/ jK \ sinni jmv*/vv" i 
/v»=/o(K)exp(-— % J /_ 2 X > (4) 

, , , p x 4 /3\3/2 u 2 ( 3u 2 \ 

/o" =— -7=(~) — exp -— 2 . (5) 

In table 1 are shown the capture rates obtained for stars at different stages of their evolution, considering different 
stellar masses M* and velocities v*, and different values for the mass of the DM particles m x . Higher capture rates 
were obtained during the MS, for stars with higher M* and lower v*, and for DM particles with lower m x . 

Once they are captured, DM particles accumulate in the center of the star, in a region with a typical radius 
r x = y / 3KT c /27iGp c m x . There, WIMPs annihilate among themselves, thus providing a new source of energy that 
contributes to the total luminosity of the star with: 

L x=fx c x m x> ( 6 ) 



where f x = 2/3 is a factor to take into account that a third of the products of DM annihilation will escape out of the 
star in the form of neutrinos. 
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FIGURE 1. Properties of stars of 1 M© that evolve in halos with DM densities from 10 5 to 10 12 GeV cm~ 3 : a) luminosity from 
the thermonuclear reactions and from DM annihilation; b) time spent by the star in the Main Sequence (MS) (from £ gmv <\%£j 
to X c < 0.001) and time until the convective core disappears; c) location of the convective and radiative zones inside the star; d) 
effective temperature of the star. The characteristics in Figures a), c) and d) are plotted at such an age that all the stars are already 
in energy equilibrium, i.e.: in the beginning of the MS for p x < 3 x 10 9 GeV cm~ 3 , and in the stationary states powered only by 
the energy from DM annihilation for higher DM densities. 



STELLAR EVOLUTION SCENARIOS WITHIN DENSE DARK MATTER HALOS 

The way DM influences the star depends on the amount of energy released by DM annihilation. In this section, 
this quantity depends mainly on the DM density on the host halo p x , given that the other parameters are fixed 
to the following fiducial values: v* = 220 km s _1 , \~ x — 270 km s~ ! , (J%,sd = 10~ 38 cm 2 , o Xi si = 10~ 44 cm 2 , and 
= 100 GeV. 

For "low" DM densities (p% < 3 x 10 9 GeV cm~ 3 for a star of 1 M Q ), the energy from DM annihilation is a 
complementary source of energy for the star (see Figure l.a)). Their main impacts on stellar evolution are: i) DM 
burning contributes, together with the thermonuclear energy, to compensate the gravitational collapse of the protostar, 
which will reach the hydrostatic equilibrium at a lower central temperature than in the classical scenario, leading to a 
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FIGURE 2. Stationary states reached by stars with masses from 0.7 to 3 M e when the energy from DM annihilation compensates 
the gravitational energy during the collapse. These equilibrium positions, where stars will remain for an indefinite time, are plotted 
for different dark matter halo densities, indicated in units of log(p^/GeV cm -3 ) at the side of each line. The grey lines are the 
classical evolutionary paths. 



lower rate of hydrogen burning and, therefore, an extension of the lifetime of the star in the MS (see Figure l.b)); ii) as 
the energy from DM annihilation is produced in a very concentrated region (r x ~ 0.03 the gradient of temperature 
in the core is very high, requiring the convective core to remain for a longer time than in the classical scenario in order 
to evacuate the extra input of energy in a more efficient way (see Figure l.b)). 

For high DM densities, the energy from the annihilation of captured DM particles compensates the gravitational 
collapse of the protostar when the central temperature is still below the threshold for thermonuclear reactions. In this 
new scenario, the star remains in an indefinite state of equilibrium with DM burning as its only source of energy. In 
Figure 2, we show the equilibrium positions in the Hertzsprung-Russell diagram reached by stars with masses between 
0.7 and 3 M that evolved in halos with different DM densities. These results show that stars "freeze" during their 
pre-MS track sooner than in our previous study [23], which was carried out with a more basic capture rate calculation. 
Regarding the energy transport mechanisms inside the star, convection is the only one capable of transporting the huge 
amount of energy through the whole star if it evolves in halos with very high DM densities. In the case of a star of 1 
M , the star is fully convective for p x > 10 10 GeV cm~ 3 . For DM densities between 10 9 and 10 10 GeV cm~ 3 , the star 
will remain with a convective core for the rest of its life (see Figure l.b) and c)). One of the direct consequences of the 
changes in the extension of the convective zones is a reduction on the effective temperature of the star T e /f. For a one 
solar mass star, the drop on T e ff is of the order of ~ 1000 K (see Figure l.c)). 



CONCLUSIONS 

The search for the unusual stars described in this manuscript appears as a promising alternative method to test and 
constrain the properties of DM candidates. In the higher horizontal axe of Figure 1, are shown the expected distances 
toward the center of our galaxy where stars with such unusual properties may be found, following the DM density 
profile of Bertone & Merrit [27]. These distances are shown to provide an indication of where to find these stars 
within the inner parsec of our galaxy. Many parameters with high uncertainties influence this prediction: the velocity 
distribution of the DM particles, their scattering cross sections with baryons and the DM density profile, among others. 



The orbit followed by these stars also plays an important role: the changes on the stellar velocity with respect to the 
halo can boost or diminish the capture rate (see table 1 and reference [22]), varying the contribution of DM burning to 
the total luminosity. We expect that with the future improvements on the observation of stars near the galactic center 
we will be able to give constraints to the DM particle models, providing an alternative approach to reveal the nature of 
Dark Matter. 
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